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Foreword 


Anyone  who  has  had  to  take  upper-division  courses  in  the  hard  sciences,  engineering,  or  applied 
mathematics  has  probably  encountered  convolution.  The  typical  methods  used  to  solve  a  convolution 
integral — the  graphical  method,  the  direct  integration  method,  or  the  use  of  Fourier  or  Laplace 
transforms — can  be  very  difficult  to  actually  set  up  and  solve,  even  when  worked  examples  are  available 
as  guides.  This  report  shows  a  method  of  solving  the  convolution  integral  that  should  make  it  easy  for 
anyone,  including  anyone  using  numerical  methods,  to  visualize  the  convolution  and  then  determine  the 
integrals  needed  to  solve  the  problem. 
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Theoretical  Development  of  a  Method  of  Evaluating  Convolutions 


The  convolution  h(x)  of  two  real- valued  functions  f(x)  and  g(x)  is  indicated  by  the  shorthand  notation 
h(x)=f(x)*g(x)  and  is  defined  by  the  integral: 

/'OO  /'OO 

h(x)  =  I  f  (a)g(x  —  a)da  =  I  f  (x  —  a)g(a)da! 

*•' — CO  J  —  CO 

Anyone  who  has  had  to  take  upper-division  courses  in  the  hard  sciences,  engineering  or  applied 
mathematics  has  probably  encountered  convolution.  The  typical  methods  used  to  solve  a  convolution 
integral — the  graphical  method,  the  direct  integration  method  or  the  use  of  Fourier  or  Laplace 
transforms — can  be  very  difficult  to  actually  set-up  and  solve,  even  when  worked  examples  are  available 
as  guides.  This  report  shows  a  method  of  solving  the  convolution  integral  that  should  make  it  easy  for 
anyone,  including  anyone  using  numerical  methods,  to  visualize  the  convolution  and  then  determine  the 
integrals  needed  to  solve  the  problem. 

Assume  we  want  to  convolve  two  functions:  f(x)  which  exists  only  on  the  closed  interval  [ai,bi];  and  g(x) 
which  exists  only  on  the  closed  interval  [a2,b2l  -  This  method  starts  by  placing  the  two  functions  f  and  g  in 
a  three-dimensional  rectangular  Cartesian  coordinate  system.  Let  Z|(x)=f(\)  and  z2(y)=g(y).  Since  the 
convolution  process  is  commutative,  it  does  not  matter  in  this  method  if  f  is  a  function  of  x  and  g  is  a 
function  of  y  or  vice-versa.  Let  z(x,y)=zj(x)z2(y)  or  z(x,y)=f(x)g(y).  So,  z(x,y)  will  only  exist  in  the 
closed  rectangular  domain  [x,y]=[a!<x<bi,a2<y<b2],  as  seen  in  Figure  1. 


Figure  1:  Domain  [ai<x<bi,a2<y<b2]  of  the  Function  z(x,y)=f(x)g(y) 


1  Linear  Systems,  Fourier  Transforms  and  Optics,  Jack  D.  Gaskill,  John  Wiley  &  Sons,  1978, 
Chapter  6,  Section  6-1,  Page  150 
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Now  translate  the  function  z(x,y)  in  the  y-direction  by  the  amount  x.  In  other  words,  move  the  rectangle 
shown  in  Figure  1  up  in  the  positive  y  direction  so  that  the  bottom  boundary  of  the  domain  is  a  distance  a2 
above  the  line  y=x.  This  changes  the  function  z(x,y)  to  the  function  z(x,y-x),  as  seen  in  Figure  2. 


Figure  2:  Domain  [ai<x<bi,a2<y-x<b2]=[ai<x<bi,x+a2<y<x+b2]  of  the  Function  z(x,y-x)=f(x)g(y-x) 

Let  a  plane  perpendicular  to  the  xy-plane  and  the  yz-plane  and  parallel  to  the  xz-plane  intersect  the 
function  z(x,y-x)  shown  in  Figure  2  at  some  point  y'e  |ai+a2,bi+b2|.  The  cross-sectional  area  formed  by 
this  plane  and  the  function  z(x,y-x)  is: 

A(y')  =  f  f(x)g(y'  -  x)dx 

Jxi 

where:  if  ai+a2<y'<ai+b2  (Area  1),  then  Xi=ai  and  x2=y'-a2 

if  ai+b2<y'<a2+bi  (Area  2),  then  X]=y'-b2  and  x2=y'-a2 
if  a2+b!<  y'<bj+b2  (Area  3),  then  X!=y'-b2  and  x2=bi 
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But,  since  the  function  z(x,y-x)  does  not  exist  outside  of  the  closed  interval  | a i , b i  ] ,  we  can  extend  the 
limits  of  integration  to  Xi=-oo  and  to  x2=+co  without  changing  the  value  of  the  integral.  This  now  gives 
the  following  equation  for  the  cross-sectional  area: 


AGO  = 


"O 

J  -  ( 


f(x)g(y' 


—  x)dx 


which  is  exactly  the  definition  of  the  convolution. 


The  following  Examples  show  how  to  apply  this  method. 


Example  1:  Convolution  of  a  Right  Triangle  Function  &  a  Rectangle  Function2 

Suppose  we  want  to  convolve  two  functions:  the  first  function  is  a  right  triangle  with  endpoints  (0,0), 
(3,0)  and  (3,2);  the  second  function  is  a  rectangle  with  endpoints  (-1,0),  (3,0),  (3,1)  and  (-1,1).  Let  the 
first  function  be  f(x)  such  that: 


f(x) 


Let  the  second  function  by  g(y)  such  that: 


0  <  x  <  3 
x  <  0,  x  >  3 


g(y) 


-1  <y<  3 
y  <  — l,y  >  3 


The  two  functions  are  shown  in  Figure  3.  Their  domain  in  the  xy-plane  is  shown  in  Figure  4,  and  their 
translated  domain  with  respect  to  the  line  y=x  is  shown  in  Figure  5 


Figure  3:  Two  Functions  to  Convolve  in  Example  1:  Right  Triangle  f(x)  and  Rectangle  g(y) 


2 

This  example  is  derived  from  the  example  shown  in  Linear  Systems,  Fourier  Transforms  and 
Optics,  Jack  D.  Gaskill,  John  Wiley  &  Sons,  1978,  Chapter  6,  Section  6-1,  Pages  154-155 
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1  1 

1  2  ■ 
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Figure  4:  Domain  |0<x<3,-l<y<3]  of  the  Function  z(x,y)=f(x)g(y) 


In  Figure  5,  there  are  three  different  sets  of  limits  of  integration  of  the  integral:  in  Area  1,  where  -l<y<2, 
the  value  of  x  can  vary  from  0  to  y+1;  in  Area  2,  where  2<y<3,  the  value  of  x  can  vary  from  0  to  3;  and  in 
Area  3,  where  3<y<6,  the  value  of  x  can  vary  from  y-3  to  3.  For  all  y<-l  and  y>6,  the  integral  is  equal  to 
zero  since  the  function  z  does  not  exist  there.  In  the  closed  y-interval  [-1,6],  we  can  evaluate  the  integral 
in  the  three  intervals  associated  with  Areas  1,  2,  and  3  to  get  a  piecewise  solution.  Plug  f(x)  into  the 
equation  for  A(y').  For  the  function  g(y),  substitute  y-x  for  y  throughout  the  function  and  then  plug  this 
function  g(y-x)  into  the  equation  for  A(y').  Since  g(y)  is  constant  in  the  interval  [-1,3],  it  will  not  change 
when  translated  from  y  to  y-x. 
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Figure  5:  Domain  |l<x<3,-l<y-x<3]=|l<x<3,x-l<y<x+3]  of  the  Function  z(x,y-x)=f(x)g(y-x) 

In  the  y-interval  associated  with  Area  1,  -l<y<2,  we  can  evaluate  the  integral  as  follows: 

/'OO 

hi(y)  =  f  (x)g(y  —  x)dx 
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(y  +  i)2 


fy  /2  \  2  fy  x 

hl(y)  =  J0  \3  /  (1)dx  =  3  J0  Xdx  =  j 

In  the  y-interval  associated  with  Area  2,  2<y<3,  we  can  evaluate  the  integral  as  follows: 

n  OO 

h2(y)  =  I  f  (x)g(y  —  x)dx 

^  —  00 

hM=l  (ix)a)dx=Uxdx=; 


=  3 


In  the  y-interval  associated  with  Area  3,  3<y<6,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h3(y)  =  I  f  (x)g(y  —  x)dx 

—  GO 


r3  (2  \  2  r3 

h3(y)  =  l-x)(l)dx  =  -  xdx  = 

^y-3  y  ^y-3 


=  3 


(y-3)2 


y— 3 


Finally,  we  can  specify  the  entire  solution  of  the  convolution  asr 


h(x)  = 


0, 


(x+1)3 


3  - 


(x-3)2 


0, 


x  <  -1 

— 1  <  x  <  2 

2  <  x  <  3 

3  <  x  <  6 
x  >  6 


Since  convolution  is  commutative,  we  will  get  the  same  answer  if  we  switch  the  functions  around  and  let 
the  first  function  be  the  rectangle  and  the  second  function  be  the  right  triangle  such  that: 


f«  =  & 


— 1  <  x  <  3 
x  <  — 1,  x  >  3 


and 


(2y 


g(y) 


_hf.  0<y<3 


0, 


y  <  0,y  >  3 

For  this  switch,  the  domain  in  the  xy-plane  is  shown  in  Figure  6,  and  the  translated  domain  with  respect  to 
the  line  y=x  is  shown  in  Figure  7. 


3  Linear  Systems,  Fourier  Transforms  and  Optics,  Jack  D.  Gaskill,  John  Wiley  &  Sons,  1978, 
Chapter  6,  Section  6-1,  Pages  154-155 
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Figure  6:  Domain  [-l<x<3,0<y<3]  of  the  Switched  Function  z(x,y)=f(x)g(y) 


In  Figure  7,  there  are  three  different  sets  of  limits  of  integration  of  the  integral:  in  Area  1,  where  -l<y<2, 
the  value  of  x  can  vary  from  0  to  y;  in  Area  2,  where  2<y<3,  the  value  of  x  can  vary  from  y-3  to  y;  and  in 
Area  3,  where  3<y<6,  the  value  of  x  can  vary  from  y-3  to  3.  For  all  y<-l  and  y>6,  the  integral  is  equal  to 
zero  since  the  function  z  does  not  exist  there.  In  the  closed  y-interval  [-1,6],  we  can  evaluate  the  integral 
in  the  three  intervals  associated  with  Areas  1,  2,  and  3  to  get  a  piecewise  solution.  Plug  f(x)  into  the 
equation  for  A(y').  For  the  function  g(y),  substitute  y-x  for  y  throughout  the  function  and  then  plug  this 
function  g(y-x)  into  the  equation  for  A(y'). 


In  the  y-interval  associated  with  Area  1,  -l<y<2,  we  can  evaluate  the  integral  as  follows: 

n  OO 

hi(y)  =  I  f (x)g(y  —  x)dx 

J  —  oo 


hi(y)  =  J  (l)  (y  —  x)) d X  =  (y - x) dx  =  ^ (xy - y) 


-1 


hi(y)  = 


2(y 


+  y  +  o  =  o(y2  +  2y  +  D  = 


(y  +  1)2 


3\2  2  J  3  "  3 

In  the  y-interval  associated  with  Area  2,  2<y<3,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h2(y)  =  I  f  (x)g(y  —  x)dx 

J  —  oo 

(•y  n  \  2  fy  2  x2 

h2(y)  =  J  (l)^-(y-x)Jdx  =  -J  (y  —  x)  dx  =  -  (xy  —  — 


h2(y)  =  3 


T2  y 


y )  -  ( y(y  -  3) 


y-3 

(y-3)2 


2  / y2 


6y  9 
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h2(y) 


Figure  7:  Domain  [-l<x<3,0<y-x<3]=|l<x<3,x<y<x+3]  of  the  Switched  Function  z(x,y-x)=f(x)g(y-x) 


In  the  y-interval  associated  with  Area  3,  3<y<6,  we  can  evaluate  the  integral  as  follows: 


h3 


f  (x)g(y  -  x)dx 
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3 


h3(y) 

h3(y) 

h3(y) 


|_a)(|(y-x))dx  =  ||_3(y-x) 

§[(3y-?)-(y(y-3)-^-) 


2  x 
dx  =  3  (xy  -  T> 


,  y: 2  ,  6y 

3y_y+y_ 


^(3y-?-v  +  o )  =  ^(3y-y2)  =  2y-  V 


By  completing  the  square,  we  can  change  this  equation  as  follows: 

h3(y)  =  -^(y2  -  6y)  =  -i(y2  -  6y  +  9  -  9)  =  - 1 ((y -  3)2  -  9)  =  3  - 


Finally,  we  can  again  specify  the  entire  solution  of  the  convolution  exactly  as  before: 

f  0,  x  <  —  1 


h(x)  = 


(x  +  l)3 

3  ' 

3, 


— 1  <  x  <  2 
2  <  x  <  3 


(x  -  3)2 

3  —  - — ,  3  <  x  <  6 

3 

y  0,  x  >  6 

This  is  exactly  the  same  as  the  previous  solution.  The  graph  of  this  solution  is  shown  in  Figure  8. 


Figure  8:  Results  of  the  Convolution  of  the  Functions  in  Example  1 


Example  2:  Convolution  of  a  Right  Triangle  Function  with  Itself 

Suppose  we  want  to  convolve  the  following  probability  density  function  (PDF)  with  itself: 

1  -  0  <  x  <  2 

2 

0,  x  <  0,  x  >  2 

This  PDF  is  a  right  triangle  with  endpoints  (0,0),  (2,0)  and  (0,1)  and  is  shown  in  Figure  9.  To  convolve 
this  PDF  with  itself,  the  second  function  will  then  be  g(y)  such  that: 

„>.M' 

(  0, 
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0  <  y  <  2 
y  <  0,y  >  2 


The  domain  of  these  two  PDF  s  in  the  xy-plane  is  shown  in  Figure  10,  and  their  translated  domain  with 
respect  to  the  line  y=x  is  shown  in  Figure  1 1 . 

In  Figure  1 1 ,  there  are  two  different  sets  of  limits  of  integration  of  the  integral:  in  Area  1 ,  where  0<y<2, 
the  value  of  x  can  vary  from  0  to  y;  in  Area  2,  where  2<y<4,  the  value  of  x  can  vary  from  y-2  to  2.  For  all 
y<0  and  y>4,  the  integral  is  equal  to  zero  since  the  function  z  does  not  exist  there.  In  the  closed  y-interval 
[0,4],  we  can  evaluate  the  integral  in  the  two  intervals  associated  with  Areas  1  and  2  to  get  a  piecewise 
solution.  Plug  f(x)  into  the  equation  for  A(y').  For  the  function  g(y),  substitute  y-x  for  y  throughout  the 
function  and  then  plug  this  function  g(y-x)  into  the  equation  for  A(y'). 


Figure  9:  The  Right  Triangle  Function  to  Convolve  with  Itself  in  Example  2 

In  the  y-interval  associated  with  Area  1,  0<y<2,  we  can  evaluate  the  integral  as  follows: 


hi(y) 

hi(y) 

hi(y) 


"O 

J  — 


f  (x)g(y  —  x)dx 


y-x 

2 


y 

4 


2y  +  - 


p  r  y+i 

y  +  x)  dx  =  -  I  (4  —  2y  +  xy  —  x2)  dx 
4  Jn 


77  (24  —  12y  +  y2) 
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Figure  10:  Domain  |0<x<2,0<y<2]  of  the  Function  z(x,y)=f(x)g(y) 
In  the  y- interval  associated  with  Area  2,  2<y<4,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h2(y)  =  f  (x)g(y  —  x)dx 


x\ 

/  y-x\ 

1/ 

(1 

(l 

1  dx  =  -  ( 

V 

22 

V  2  > 

4  \ 

f 

Jy-2 


1  /  g 

h2  (y)  =  — f8  —  4y  +  2y  —  —  —  (4(y  —  2)  —  2y(y  -  2)  + 


y— 2 

y(y-2)2  (y-2)3 


h2(y)=^(y-2y-(4y-8-2y2+4y  +  y 


2  3 

3  -  4y2  +  4y  y3  -  6y2  +  12y  -  8 


h2(y)  =t 


2  3 

1  / 32  -  60y  +  48  +  12y2  -  3y3  +  12y2  -  12y  +  2y3  -  12y2  +  2y3  -  12y2  +  24y  -  16 


h2(y)  = 


64  -  48y  +  12y2  -  y3 
24 
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Figure  11:  Domain  |0<x<2,0<y-x<2]=|0<x<2,x<y<x+2]  of  the  Function  z(x,y-x)=f(x)g(y-x) 

Finally,  we  can  specify  the  entire  solution  of  the  convolution  as: 

f  0,  x  <  0 


h(x)  = 


x 


24(24-  12x  +  x2), 
64  -  48x  +  12x2  -  x3 


24 


0, 


0  <  x  <  2 

2  <  x  <  4 
x  >  4 


The  graph  of  this  solution  is  shown  in  Figure  12. 
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Figure  12:  Results  of  the  Convolution  of  the  Function  in  Example  2  with  Itself 


Another  advantage  of  this  method  of  solving  convolutions  by  direct  integration  is  that  it  can  allow  people 
to  solve  convolutions  of  functions  that  have  variable  shapes,  provided  that  the  functions  being  convolved 
are  independent  of  each  other.4  The  following  example  shows  how  to  do  this. 


Example  3:  Convolution  of  a  Right  Triangle  Function  &  a  Varying  Isosceles  Triangle 
Function 

Suppose  we  want  to  find  the  convolution  of  two  probability  density  functions  (PDFs).  Let  the  first  PDF 
be: 

3  —  x 

- ,  1  <  x  <  3 

2 

0,  x  <  1,  x  >  3 

This  PDF  is  a  right  triangle  with  vertices  at  (1,0),  (1,1)  and  (3,0).  Let  the  second  PDF,  g(x),  be  an 
isosceles  triangle  centered  (and  having  its  mode)  at  any  point  on  the  positive  x-axis.  The  base  of  this 
triangle  extends  to  ±50%  of  the  central  value  on  the  x-axis,  and  the  probability  goes  to  zero  at  the  ends  of 
this  base.  Since  the  length  of  the  base  increases  as  the  central  value  on  the  x-axis  increases,  the  height  of 
the  triangle  must  decrease  as  the  central  value  on  the  x-axis  increases  to  keep  the  area  of  the  PDF  equal  to 


4  An  example  of  a  varying  shape  would  be  the  PDF  associated  with  a  multiplicative  estimating 
error.  Such  an  error  is  typically  stated  as  a  plus-or-minus  percent  of  a  given  value  with  the  PDF 
mode  at  the  given  value  and  the  PDF  going  to  zero  at  the  maximum  plus  or  minus  percent. 
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one.  So,  the  length  of  the  base,  b,  of  the  isosceles  triangle  PDF  centered  at  a  point  xc  on  the  x-axis  will  be 
b=2(xc/2)=xc.  Using  the  standard  equation  for  the  area  of  a  triangle,  we  can  then  calculate  the  height  of 
the  isosceles  triangle  PDF : 


To  find  the  convolution  of  these  two  PDFs,  keep  the  first  function,  the  right -triangle  PDF,  as  a  function  of 
x,  f(x).  Make  the  second  function,  the  isosceles  triangle  PDF  a  function  of  y,  g(y).  For  the  second  PFD, 
the  isosceles  triangle  is  still  centered  at  any  point  x  on  the  positive  x-axis,  but  this  triangle  is  now  in  a 
plane  parallel  to  the  yz-plane  and  perpendicular  to  the  xz-plane  with  endpoints  at  (x,x/2,0),  (x,-x/2,0)  and 
(x,0,2/x).  For  the  side  of  the  triangle  in  the  positive  y  region,  the  slope  of  the  line  segment  formed  by  this 
side  is  m+=(2/x-0)/(0-x/2)=(-4/x2),  and,  for  the  side  of  the  triangle  in  the  negative  y  region,  the  slope  of 
the  line  segment  formed  by  this  side  is  m_=(2/x-0)/(0-(-x/2))=(4/x2).  So,  the  equation  of  the  line  segment 
in  the  positive  y  region  is  z-0=(-4/x2)(y-x/2)=-2(2y-x)/x2,  and  the  equation  of  the  line  segment  in  the 
negative  y  region  is  z-0=(4/x2)(y-(-x/2))=2(2y+x)/x2. 


The  two  functions  are  shown  in  Figure  13.  The  domain  of  these  two  functions  in  the  xy-plane  is  shown  in 
Figure  14,  and  Figure  15  shows  what  the  function  z(x,y)=f(x)g(y)  looks  like.  The  translated  domain  of  the 
two  functions  with  respect  to  the  line  y=x  is  shown  in  Figure  16. 


Figure  13:  The  Two  Functions  to  Convolve  in  Example  3 
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Figure  14:  Domain  |l<x<3,-x/2<y-x<x/21=|l<x<3,x/2<y<3x/2]  of  the  Function  z(x,y-x)=f(x)g(y-x) 


Since  g(y)  is  a  piece-wise  continuous  function  described  by  one  function  when  y>0  and  a  different 
function  when  y<0,  then  g(y-x)  is  also  a  piece-wise  continuous  function  described  by  one  function — call 
it  g(y-x)y>x — when  y>x  and  a  different  function — call  it  g(y-x)y<x — when  y<x.  Also,  as  shown  in  Figure 
16,  there  are  four  different  sets  of  limits  of  integration  of  the  integrals:  in  Area  1,  where  0.5<y<l,  the 
value  of  x  can  vary  from  1  to  2y  for  the  integrand  f(x)g(y-x)y<x;  in  Area  2,  where  l<y<1.5,  the  value  of  x 
can  vary  from  1  to  y  for  the  integrand  f(x)g(y-x)y>x  (in  Subarea  2,1)  and  from  y  to  2y  for  the  integrand 
f(x)g(y-x)y<x  (in  Subarea  2,2);  in  Area  3,  where  1.5<y<3,  the  value  of  x  can  vary  from  2y/3  to  y  for  the 
integrand  f(x)g(y-x)y>x  (in  Subarea  3,1)  and  from  y  to  3  for  the  integrand  f(x)g(y-x)y<x  (in  Subarea  3,2); 
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and  in  Area  4,  where  3<y<4.5,  the  value  of  x  can  vary  from  2y/3  to  3  for  the  integrand  f(x)g(y-x)y>x.  For 
all  y<0.5  and  y>4.5,  the  integral  is  equal  to  zero  since  the  function  z  does  not  exist  there.  In  the  closed  y- 
interval  [0.5, 4. 5],  we  can  evaluate  the  integral  in  the  four  intervals  associated  with  Areas  1,  2,  3  and  4  to 
get  a  piecewise  solution.  Plug  f(x)  into  the  equation  for  A(y').  For  the  function  g(y),  substitute  y-x  for  y 
throughout  the  function  and  then  plug  this  function  g(y-x)  into  the  equation  for  A(y'). 


1.8-2 
1.6-1.8 
1.4-1.6 
1.2-1.4 
I  1-1.2 
I  0.8-1 
l  0.6-0. 8 
l  0.4-0.6 
I  0. 2-0.4 
l  0-0.2 


Figure  15:  The  Surface  of  the  Product  of  the  Two  Functions  z(x,y)=f(x)g(y)  in  Example  3 

In  the  y- interval  associated  with  Area  1,  0.5<y<l,  we  can  evaluate  the  integral  as  follows: 

/'OO 

hi(y)  =  I  f (x)g(y  —  x)dx 

—  GO 


2y  /3  —  x\  / 2(2(y-x)  +x) 


dx 


-f 


(3  -  x)  (2y  —  x) 


dx 


2y 


hi(y)  =  J  -  3  ^2y  +  l)dx=  (~y_(3  +  2y)  lnW  +  x) 

hi(y)  =  (-3  -  (3  +  2y)  ln(2y)  +  2y)  -  (-6  -  0  +  1)  =  8y  -  (3  +  2y)  ln(2y)  -  4 
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Figure  16:  Domain  |l<x<3,-x/2<y-x<x/21=|l<x<3,x/2<y<3x/2]  of  the  Function  z(x,y-x)=f(x)g(y-x) 


In  the  y-interval  associated  with  Area  2,  l<y<1.5,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h2(y)  =  |  f  (x)g(y  —  x)dx 

— QO 


h2(y)  =  J 


3  —  x\  /  2  (2  (y  —  x)  —  x) 


dx  + 


2y/3-x\  /2(2(y  —  x)  +  x) 


/3  -  x\ 


■ 


dx 


(x-  3)(2y-  3x) 


dx  + 


I 


2y 


(3  —  x)(2y  —  x) 


dx 
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■f^X 


.  ^  fy  (  6y  9  +  2y  \  f  zy  /6y  3  +  2y  \ 

h2<y)=l(~^+— — 3)dx+Jy  b — — +1)dx 

/6y  \  y  /  6y  \  2y 

h2(y)  =  +  (9  +  2y)  ln(x)  -  3xJ  ^  +  (y  —  -  (3  +  2y)  ln(x)  +  xj 

h2(y)  =  (6  +  (9  +  2y)  ln(y)  -  3y)  -  (6y  -  3)  +  (-3  -  (3  +  2y)  ln(2y)  +  2y) 

-  (-6  -  (3  +  2y)  ln(y)  +  y) 

h2(y)  =  12  -  8y  +  4(3  +  y)  ln(y)  -  (3  +  2y)  ln(2y) 

In  the  y-interval  associated  with  Area  3,  1.5<y<3,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h3(y)=  f  (x)g(y  —  x)dx 

J  —  oo 

[y  (3-x\(  2(2(y-x)-x)'\  (3  {3  -  x\  (2(2(y  -  x)  +  x)\ 

h3(y)=jfi— J( — ^ — )dx+jy  t— J( — ^ — r 

,  N  fy  (x  -  3)(2y  -  3x)  f 3  (3  -  x)(2y  -  x) 

h3(y)  =  k - ^ - dx + i - ^ 


+  1  dx 


3-x\  / 2  (2  (y  —  x)  +x) 


(3  ~  X) 
kv  \~1~J 


(x  3)(2y  3x) 


(3  —  x)(2y  —  x) 


,  ,  ,  fy(  6y  9  +  2y  \  f 3  (Gy 

y)  =  Ja V  x*  3  x - 3)dx  +  l  fe- 


6y  3  +  2y 


+  dx 


h3(y)  =  (Y  +  (9  +  2y)ln(x)_3x)|2y  +  (“V_(3  +  2y)ln(x)  +  X)l 

“3”  y 

h3(y)  =  (6  +  (9  +  2y)  ln(y)  -  3y)  -  (9  +  (9  +  2y)  In  (y)  -  2y)  +  (-2y  -  (3  +  2y)  ln(3)  +  3) 
-  (-6  -  (3  +  2y)  ln(y)  +  y) 

h3(y)  =  4(3  +  y)  ln(y)  -  (9  +  2y)  In  (|)  -  (3  +  2y)  ln(3)  +  6  -  4y 

In  the  y-interval  associated  with  Area  4,  3<y<4.5,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h4(y)  =  I  f  (x)g(y  —  x)dx 

*'  — OO 

b* w  =  k  (• 3~r)  (“  —  f  -)  dx  =  (?  +  (9  +  2y)  ln(x)  “ 3x)  la 

3  V  '  3 

h4(y)  =  (2y  +  (9  +  2y)  ln(3)  -  9)  -  (9  +  (9  +  2y)  In  (y)  -  2y) 
h4(y)  =  (9  +  2y)  In  (^j  +  4y  -  18 

Finally,  we  can  specify  the  entire  solution  of  the  convolution  as: 
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h(x)  = 


0, 

8x  —  (3  +  2x)  ln(2x)  —  4, 
12  —  8x  +  4(3  +  x)  ln(x)  —  (3  +  2x)  ln(2x) , 
2xn 


4(3  +  x)  ln(x)  —  (9  +  2x)  In  j  —  (3  +  2x)  ln(3)  +  6  —  4x, 

9 


(9  +  2x)  In  J  +  4x 


18, 

0, 


The  graph  of  this  solution  is  shown  in  Figure  17. 


x  <  0.5 
0.5  <  x  <  1 
1  <x<  1.5 

1.5  <  x  <  3 

3  <  x  <  4.5 
x  >  4.5 


Figure  17:  The  Results  of  the  Convolution  of  the  Functions  in  Example  3 
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Repeated  Convolution 


Given  the  convolution  h(x)=f(x)*g(x),  it  is  generally  true  that  h(x)  will  be  smoother  than  either 
f(x)  or  g(x),  and  the  convolution  of  a  large  number  of  functions  generally  yields  a  function  that  is 
much  smoother  than  any  of  the  individual  functions  involved  in  the  convolution.5  In  addition,  as 
the  number  of  functions  convolved  becomes  larger  and  larger,  the  resulting  convolution  often 
begins  to  look  more  and  more  like  a  Gaussian  function.  For  example,  consider  the  convolution: 

h(x)  =  fi  (x)  *  f2(x)  *  f3(x)  *  ...  *  fn(x) 

It  can  be  shown  that,  if  the  functions  fi(x),  f2(x),  f3(x),...,fn(x)  satisfy  certain  conditions,  then  as 
n— >oo,  h(x)  will  tend  towards  a  Gaussian  function.  This  is  a  non-rigorous  statement  of  the 
Central  Limit  Theorem.6  Example  4  shows  these  behaviors. 


Example  4:  Repeated  Convolution  of  a  Rectangle  Function  with  Itself 

Suppose  we  want  to  find  the  results  of  the  repeated  convolution  of  the  following  rectangle  PDF: 


1  <  x  <  2 
x  <  1,  x  >  2 


The  PDF  is  shown  in  Figure  18  along  with  its  corresponding  Gaussian  function,  i.e.,  a  Gaussian 
function  with  the  same  mean  and  variance  as  a  given  PDF.  Addendum  1  shows  how  to  calculate 
corresponding  normal  and  lognormal  functions. 


Figure  18:  Rectangle  Function  in  Example  4  &  Corresponding  (Dotted)  Gaussian  Function 


To  convolve  this  PDF  with  itself,  the  second  function  will  then  be  g(y)  such  that: 


1  <y<  2 
y  <  i.y  >  2 


The  domain  of  these  two  functions  in  the  xy-plane  is  shown  in  Figure  19,  and  their  translated  domain  with 
respect  to  the  line  y=x  is  shown  in  Figure  20. 


5  Linear  Systems,  Fourier  Transforms  and  Optics,  Jack  D.  Gaskill,  John  Wiley  &  Sons,  1978, 
Chapter  6,  Section  6-3,  Page  162,  Paragraph  2  and  Page  163  Paragraph  2. 

6  Ibid,  Page  164,  Paragraph  1. 
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Figure  19:  Domain  [l<x<2,l<y<2]  of  the  Function  z(x,y)=f(x)g(y) 


In  Figure  20,  there  are  two  different  sets  of  limits  of  integration  of  the  integral:  in  Area  1,  where  2<y<3, 
the  value  of  x  can  vary  from  1  to  y-1;  in  Area  2,  where  3<y<4,  the  value  of  x  can  vary  from  y-2  to  2.  For 
all  y<2  and  y>4,  the  integral  is  equal  to  zero  since  the  function  z  does  not  exist  there.  In  the  closed  y- 
interval  [2,4],  we  can  evaluate  the  integral  in  the  two  intervals  associated  with  Areas  1  and  2  to  get  a 
piecewise  solution.  Plug  f(x)  into  the  equation  for  A(y').  For  the  function  g(y),  substitute  y-x  for  y 
throughout  the  function  and  then  plug  this  function  g(y-x)  into  the  equation  for  A(y'). 

In  the  y-interval  associated  with  Area  1,  2<y<3,  we  can  evaluate  the  integral  as  follows: 

n  GO 

hi(y)  =  I  f (x)g(y  —  x)dx 

J — oo 

hi(y)  =  J  (1)  (l)dx  =  J  dx  =  x|^_1  =  y— 1  —  l  =  y  —  2 

In  the  y-interval  associated  with  Area  2,  3<y<4,  we  can  evaluate  the  integral  as  follows: 
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h2(y) 

h2(y) 


/'OO 

I  f  (x)g(y  -  x)dx 

J  —  GO 

2  2 

[  (1)  (l)d x=  [  dx  =  x|y_2  =  2  —  (y  —  2) 


4  -  y 


Figure  20:  Domain  |l<x<2,l<y-x<2]=  |l<x<2,x+l<y<x+2]  of  the  Function  z(x,y-x)=f(x)g(y-x) 

The  entire  solution  of  the  first  convolution  is: 

(x-  2, 

h(x)  =  |4-x, 

1  0, 
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2  <  x  <  3 

3  <  x  <  4 

y  <  2,y  >  4 


This  function  is  shown  in  Figure  21  along  with  its  corresponding  Gaussian  function,  which  makes 
a  fairly  good  fit  to  the  convolution  function.  If  two  functions  are  not  correlated,  i.e.,  if  they  are 
independent  functions,  then  the  mean  of  the  convolution  is  the  sum  of  the  means  of  the  two 
functions  being  convolved,  and  the  variance  of  the  convolution  is  the  sum  of  the  variances  of  the 
two  functions  being  convolved.  Addendum  2  shows  the  proofs  of  these  statements.  So,  for 
repeated  convolutions,  the  mean  and  variance  of  the  corresponding  Gaussian  function  is  just  the 
respective  sums  of  the  means  and  variances  of  the  functions  being  convolved. 


Figure  21:  Results  of  1st  Example  4  Convolution  &  Corresponding  (Dotted)  Gaussian  Function 


For  the  second  convolution,  let  this  function  h(x)  from  the  first  convolution  be  the  new  f(x)  and 
convolve  it  with  the  original  g(y).  In  this  case,  f(x)  is  a  piece-wise  continuous  function  described  by 
one  function  when  2<x<3 — call  it  fi(x) — and  a  different  function  when  3<x<4 — call  it  f2(x).  The  domain 
of  these  two  functions  in  the  xy-plane  is  shown  in  Figure  22,  and  their  translated  domain  with  respect  to 
the  line  y=x  is  shown  in  Figure  23. 


In  Figure  23,  there  are  three  different  sets  of  limits  of  integration  of  the  integral:  in  Area  1,  where  3<y<4, 
the  value  of  x  can  vary  from  2  to  y-1  for  the  integrand  fi(x)g(y-x);  in  Area  2,  where  4<y<5,  the  value  of  x 
can  vary  from  y-2  to  3  for  the  integrand  fj(x)g(y-x)  (in  Subarea  2,1)  and  from  3  to  y-1  for  the  integrand 
f2(x)g(y-x)  (in  Subarea  2,2);  and  in  Area  3,  where  5<y<6,  the  value  of  x  can  vary  from  y-2  to  4  for  the 
integrand  f2(x)g(y-x).  For  all  y<3  and  y>6,  the  integral  is  equal  to  zero  since  the  function  z  does  not  exist 
there.  In  the  closed  y-interval  [3,6],  we  can  evaluate  the  integral  in  the  three  intervals  associated  with 
Areas  1,  2  and  3  to  get  a  piecewise  solution.  Plug  f(x)  into  the  equation  for  A(y').  For  the  function  g(y), 
substitute  y-x  for  y  throughout  the  function  and  then  plug  this  function  g(y-x)  into  the  equation  for  A(y'). 


In  the  y-interval  associated  with  Area  1,  3<y<4,  we  can  evaluate  the  integral  as  follows: 
hi(y)  =  I  f  (x)g(y  —  x)dx 

j  —  00 


hi(y)  =  j 

hi(y) 


y-1 


fx2  \ 

y-1 

2 

(y-D2 


—  2(y— 1) 


-[2-4] 


y2  —  2y  +  1  n  ,  n  ,  n  y2- 6y  +  9  (y-3)2 

- - - 2y  +  2  +  2  = - - - = - - - 
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Figure  22:  Domain  [2<x<4,l<y<2]  of  the  Function  z(x,y)=f(x)g(y) 


In  the  y-interval  associated  with  Area  2,  4<y<5,  we  can  evaluate  the  integral  as  follows: 


/'OO 

h2(y)  =  I  f  (x)g(y  —  x)dx 

— oo 


(x— 2)(l)dx+  I  (4  — x)(l)dx=  — — 2x  +  4x — — 


h2(y)  =  [  (x  —  2)  (l)dx  4 

Jy—2 

r  A  r9  [(y-2)2 

h2(y)=  ^-6 - — 


2(y  —  2)  +  4(y  —  1) 


(y  - 1): 


9 

12  “2 


y2  -  4y  +  4  y2  -  2y  +  1 

h2(y)  =  9  —  18  —  - - ^ - +  2y  —  4  +  4y  —  4  —  - - -j- - 


h2(y)  =  -17  +  6y- 


2y2  -  6y  +  5  -34  +  12y  -  2y2  +  6y  -  5  2y2  -  18y  +  39 


In  the  y-interval  associated  with  Area  3,  5<y<6,  we  can  evaluate  the  integral  as  follows: 


/'OO 

h3(y)  =  I  f  (x)g(y  —  x)dx 

h3(y)  =  [  (4  —  x)  (l)dx  = 

K-2 


(4  —  x)  (l)dx  =  (  4x_y  )  =  [16-8]  -  4(y  —  2) 


(y-2)2 

2 


h3(y)  =  8-4y  +  8  + 


y2  -  4y  +  4  y2  -  12y  +  36  (y  -  6)2 
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Figure  23:  Domain  [2<x<4,l<y-x<2]=|2<x<4,x+l<y<x+2)  of  the  Function  z(x,y-x)=f(x)g(y-x) 


The  entire  solution  of  the  second  convolution  is: 

'  (y-3)2 


h(x)  = 


2y2  -  18y  +  39 

2  ’ 
(y-6)2 


3  <  x  <  4 

4  <  x  <  5 

5  <  x  <  6 


v  0,  y  <  3,y  >  6 

This  function  is  shown  in  Figure  24  along  with  its  corresponding  Gaussian  function,  which  makes 
a  good  fit  to  the  convolution  function. 
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0  1  2  3  4  5  6  7  X 

Figure  24:  Results  of  2nd  Example  4  Convolution  &  Corresponding  (Dotted)  Gaussian  Function 

For  the  third  convolution,  let  this  function  h(x)  from  the  second  convolution  be  the  new  f(x)  and 
convolve  it  with  the  original  g(y).  In  this  case,  f(x)  is  a  piece-wise  continuous  function  described  by 
one  function  when  3<x<4 — call  it  fj(x) — a  second  function  when  4<x<5 — call  it  f2(x) — and  a  third 
function  when  5<x<6 — call  it  f3(x).  The  domain  of  these  three  functions  in  the  xy-plane  is  shown  in 
Figure  25,  and  their  translated  domain  with  respect  to  the  line  y=x  is  shown  in  Figure  26. 

In  Figure  26,  there  are  four  different  sets  of  limits  of  integration  of  the  integral:  in  Area  1,  where  4<y<5, 
the  value  of  x  can  vary  from  3  to  y-1  for  the  integrand  fi(x)g(y-x);  in  Area  2,  where  5<y<6,  the  value  of  x 
can  vary  from  y-2  to  4  for  the  integrand  fi(x)g(y-x)  (in  Subarea  2,1)  and  from  4  to  y-1  for  the  integrand 
f2(x)g(y-x)  (in  Subarea  2,2);  in  Area  3,  where  6<y<7,  the  value  of  x  can  vary  from  y-2  to  5  for  the 
integrand  f2(x)g(y-x)  (in  Subarea  3,1)  and  from  5  to  y-1  for  the  integrand  f3(x)g(y-x)  (in  Subarea  3,2);  and 
in  Area  4,  where  7<y<8,  the  value  of  x  can  vary  from  y-2  to  6  for  the  integrand  f3(x)g(y-x).  For  all  y<4 
and  y>8,  the  integral  is  equal  to  zero  since  the  function  z  does  not  exist  there.  In  the  closed  y-interval 
[4,8],  we  can  evaluate  the  integral  in  the  four  intervals  associated  with  Areas  1,  2,  3  and  4  to  get  a 
piecewise  solution.  Plug  f(x)  into  the  equation  for  A(y').  For  the  function  g(y),  substitute  y-x  for  y 
throughout  the  function  and  then  plug  this  function  g(y-x)  into  the  equation  for  A(y'). 
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Figure  25:  Domain  [3<x<6,l<y<2]  of  the  Function  z(x,y)=f(x)g(y) 
In  the  y- interval  associated  with  Area  1,  4<y<5,  we  can  evaluate  the  integral  as  follows: 


n  GO 

hi(y)  =  I  f  (x)g(y  —  x)dx 

*'  — GO 

hi(y)=r  (^-)° 


1  fy_1  1  /x3 

(l)dx  =  -  (x2  —  6x  +  9)  dx  =  - 1  — —  3x2  +  9x 

2  Jo  2  V  3 


hi(y)  =\i  ^  ~  3(y  —  l)2  +  9(y  —  1)  -[9  -  27  +  27] 


1  (\y3  -  3y2  +  3y  -  1  -  9y2  +  18y  -  9  +  27y  -  27 


hi(y)  = 


hi(y)  = 


In  the  y-interval  associated  with  Area  2,  5<y<6,  we  can  evaluate  the  integral  as  follows: 


1  /y3  -  12y2  +  48y  -  64\  (y  -  4) 


/'OO 

h2(y)  =  |  f  (x)g(y  —  x)dx 

*'  — oo 

hM  =  L(^ir-)(1) 


(l)dx  + 


y_1  f  2x2  -  18x  +  39 


(l)dx 


1  /  2x3 


1  /x3  \  1  /2x3  \  y_ 

h2Cy)  =  2  (y  “  3x2  +  9x)  y_2  “  2  (“T  ~  9x2  +  39x)  ^ 

h2(y)  =  1  ([^  -  48  +  36]  -  [(y-^-ocy-a^^-a- 
[2(y  —  l)3  -  27(y  —  l)2  +  117(y  -  1)1  [128 


144  + 156 
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(192  -  y3  +  6y2  -  12y  +  8  +  9y2  -  36y  +  36  -  27y  +  54  -  2y3  +  6y2  -  6y  +  2  +  27y 
—  54y  +  27  —  117y  +  117) 

1 

-  (3y3  -  48y2  +  252y  -  436) 


y=X+Z  /jSubare 

x=y-2//'Area  3/ 


Area2/ry=x+l 


( 

i  X 

0  12  3 


Figure  26:  Domain  [3<x<6,l<y-x<2]=|3<x<6,x+l<y<x+21  of  the  Function  z(x,y-x)=f(x)g(y-x) 

In  the  y-interval  associated  with  Area  3,  6<y<7,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h3(y)  =  |  f  (x)g(y  —  x)dx 

—  GO 

y-Vfx-6)2 


9x2  +  39x 


6x2  +  36x 


h3(y)  = 


250 


225  + 195 


+ 


2(y  —  2)3  —  27(y  —  2)2  +  117(y  —  2) 


+ 


(y-l)3-18(y-l)2  +  108(y-l) 


125 


150  + 180 


h3(y)  =  -  (2y3  -  12y2  +  24y  -  16  -  27y2  +  108y  -  108  +  117y  -  234  +  y3  -  3y2  +  3y  -  1 


18y2  +  36y  -  18  +  108y  -  375) 


h3(y)  =  -  (3y3  -  60y2  +  396y  -  860) 

6 

In  the  y-interval  associated  with  Area  4,  7<y<8,  we  can  evaluate  the  integral  as  follows: 

/'OO 

h4(y)  =  I  f  (x)g(y  —  x)dx 

*'  — OO 

(y)  =  f 

Jy-2  , 

(y-2)3 


(x  —  6)2\  1/x3 

(1)dx=2  y~6x  +36x 


y-2 


h4(y)  =  -  (  [72  -  216  +  216]  - 


—  6(y  —  2)2  +  36(y  —  2) 


h4(y)  =  -  (216  -  y3  +  6  -  12y  +  8  +  18y2  +  72y  -  72  +  108y  -  216) 
6 


h4(y)  =  -  7 


1  fy3  —  24y2  +  192y-512 
6 


(y-  8)2 


3  J  6 

The  entire  solution  of  the  third  convolution  is: 

r  (y-4)3 


h(x)  = 


-  -  (3y3  -  48y2  +  252y  -  436), 
6 
1 

-  (3y3  -  60y2  +  396y  -  860), 
(y-8)3 


0, 


4  <  x  <  5 

5  <  x  <  6 

6  <  x  <  7 

7  <  x  <  8 

y  <  4,  y  >  8 


This  function  is  shown  in  Figure  27  along  with  its  corresponding  Gaussian  function,  which  makes 
an  extremely  good  fit  to  the  convolution  function. 

Because  the  rectangle  function  in  Example  4  is  a  symmetrical  function  about  the  line  x=3/2,  the 
repeated  convolutions  of  this  rectangle  function  converged  quite  rapidly — in  just  a  few 

-7 

convolutions — to  their  corresponding  Gaussian  function  .  But,  the  repeated  convolution  of 
highly  asymmetrical  functions,  as  shown  in  Example  5,  may  take  much  longer  to  converge  well 
to  a  corresponding  Gaussian  function. 


7  Linear  Systems,  Fourier  Transforms  and  Optics,  Jack  D.  Gaskill,  John  Wiley  &  Sons,  1978, 
Chapter  6,  Section  6-3,  Page  165,  Paragraph  2 
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Figure  27:  Results  of  3rd  Example  4  Convolution  &  Corresponding  (Dotted)  Gaussian  Function 


Example  5:  Repeated  Convolution  of  an  Exponential  Function  with  Itself8 

Suppose  we  want  to  find  the  results  of  the  repeated  convolution  of  the  following  exponential  PDF: 


Then,  g(y)  is: 


And,  g(y-x)  is: 


x  >  0 
x  <  0 


y>  0 
y  <  0 


g(y  -  x)  =  {' 


e-(y-x); 
0, 


y  —  x  >  0 
y  —  x  <  0 


The  domain  of  the  product  f(x)g(y)  in  the  xy-plane  is  the  entire  first  quadrant,  and  the  translated  domain 
with  respect  to  the  line  y=x  of  the  product  f(x)g(y-x)  is  the  area  in  the  first  quadrant  between  the  line  y=x 
and  the  y-axis  as  shown  in  Figure  28. 


o 

This  example  is  derived  from  the  example  shown  in  Linear  Systems,  Fourier  Transforms  and 
Optics,  Jack  D.  Gaskill,  John  Wiley  &  Sons,  1978,  Chapter  6,  Section  6-3,  Page  165 
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Figure  28:  Domain  [x>0,y-x>0]=[x>0,y>x]  of  the  Function  z(x,y-x)=f(x)g(y-x) 

The  y-interval  for  this  entire  domain  is  [0,<x>),  so  we  can  evaluate  the  integral  for  hn(y),  the  nth  convolution 
of  f(x),as  follows: 

n  GO 

hi(y)  =  |  f(x)g(y-x)dx 
j  —  00 

ry  ry  ry 

hi(y)  =  (  (e-x)  (e_(^y_x'))dx  =  I  e_x_y+x  dx  =  I  e_y  dx  =  xe-y|Q  =  ye_y 

J  0  ”'0  ”'0 

If  we  repeat  the  convolution  by  always  letting  f(x)=hn(x),  for  all  n,  and  using  g(y)=e'y,  we  will  get: 

nCO 

h2(y)  =  |  f  (x)g(y  —  x)dx 

fy  fy  fy  x2 

h2(y)  =  I  (xe_x)  (e_(^y_x-))dx  =  I  xe_x_y+xdx=  I  xe_y  dx  =  —  e_y 

j  0  ”'0  ”'0  2 
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yv2  x 

e-y  dx  =  —  e_y 
6 


n  GO 

h3(y)  =  I  f  (x)g(y  —  x)dx 

J  —  oo 

nx2  \  ryx2  ryxz 

_e-xj(e-(y-x))dx=  j  _e-x-y+xdx=  j  _ 

/'OO 

h4(y)  =  I  f  (x)g(y  —  x)dx 

J  —  co 

nx3  \  fyx3  ryx3 

—  e~xj(e~<-y~x^)dx=J  —  e~x~y+x  dx  =  J  —e~ydx  = 


24 


n-y 


y  3 
y 


y  y4 

=  —  e“y 
o  24 


From  these  first  four  repeated  convolutions,  the  resulting  convolution  seems  to  follow  the 


following  pattern,  such  that  the  nth  convolution  will  be: 
yn 

hn(y)=^!e_y 


where  n=0,  1,  2,...  (i.e.,  n  can  be  all  non-negative  integers).  If  we  plug  in  0,  1,  2,  3  and  4  for  n  in 
this  equation  hn(y),  we  get  the  same  results  as  we  calculated  above.  To  show  this  is  true  in 
general,  we  can  complete  a  simple  proof-by-induction.  Assume  this  equation  h„(y)  is  true  and 
then  show  it  will  still  be  true  for  hn+i(y): 


hn+i(y) 

h„+i(y) 

hn+i(y) 


/'OO 

I  hn  (x)g(y  -  x)dx 

j  — CO 


e-x-y+x  dx 


=  [y*-, 
J0  n! 


e  y  dx 


7n+l 


(n+l)! 


3-y 


xn+1 

n!  (n  +  1) 


e“y 


y 

o 


So  the  equation  for  hn(y)  is  true  for  all  n.  The  results  of  the  first  four  repeated  convolutions  of 
the  function  in  this  example  are  shown  in  Figure  29  and  Figure  30:  Figure  29  shows  the 
corresponding  normal  functions;  and  Figure  30  shows  the  corresponding  lognormal  functions.  Also, 
Figure  31,  Figure  32  and  Figure  33  show  the  results  after  10,  20  and  30  repeated  convolutions, 
respectively.  For  the  asymmetrical  function  used  in  this  example,  the  convolutions  have  to  be 
repeated  many  times  before  they  truly  come  close  to  corresponding  normal  curves.  Actually,  the 
corresponding  lognormal  curves  seem  to  be  better  approximations,  but,  eventually,  after  a  few 
dozen  repetitions  of  the  convolution,  both  the  convolution  and  its  corresponding  lognormal  curve 
start  getting  fairly  close  to  the  corresponding  normal  curve. 
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Figure  29:  Example  5  Function,  4  Repeated  Convolutions  &  Corresponding  Normal  Functions 


Figure  30:  Example  5  Function,  4  Repeated  Convolutions  &  Corresponding  Fognormal  Functions 


38 


20th_Convolution 
—  —  CorrespondNormal 
. CorrespondLognormal 


Figure  32:  20th  Repeated  Convolution  &  Corresponding  Normal  &  Lognormal  Functions 


39 


Addendum  1:  Determining  the  Corresponding  Normal  &  Lognormal  PDF 


Determining  the  Corresponding  Normal  PDF 

To  find  the  corresponding  normal  PDF  that  has  the  same  mean,  p,  and  variance,  a2,  of  a  given  PDF,  we 
can  start  by  writing  a  normal  equation  in  the  general  form: 

f(x)  =  ae-b(-x-c')2  (1) 

where  c  and  a  are  the  abscissa  and  ordinate,  respectively,  of  the  mode  of  this  corresponding  normal  PDF. 
For  a  normal  equation,  the  mode,  median  and  mean  are  identically  the  same,  so  we  can  use  the  above 
equation  to  get: 


c  =  p 

To  find  the  other  values,  start  by  calculating  the  area  under  equation  ( 1 )  by  evaluating  the  integral 


A 


-£ 


ae  b(x  c)2 


using  the  transformation  m=x-c,  with  dm=dx,  and  where  x=-co— »m=-co  and  x=co— »m=oo.  The 
transformed  integral  is : 


A 


-■£ 


dm 


The  easiest  way  to  evaluate  A  is  by  using  the  standard  trick  from  upper-division  Calculus  in  which  we 
write  A  as  a  two-dimensional  integral  in  the  Cartesian  mn-plane  and  then  introduce  polar  coordinates: 


A2 


e 

O 

"  2ll  r  CO 


(/’OU  \  /  p  OO  \  p  (X)  p 

a  e-bm2dmj|a  e_bn2dnj  =  a2  e-b(m2+n2)  dmdn 

j  —  00  /  \  J  —  00  /  — GO  *■'  —  GO 


fzlt  f”  ,  ,  f°°  ,  ,  na2  f°°  ,  , 

A2  =  a2  I  I  e  br  rdrdB  =  2na2  I  e  br  rdr  =  — —  I  e  br  (— 2br)dr 

J  0  ”'0  -'o  b  J0 


n  Tia  ,  2 

A2  =  — —  e“br 


na2  na2 

—  (0-1)  =  — 
b  b 


A  =  Nb 

If  our  corresponding  normal  curve  represents  a  normalized  PDF,  then  A=l,  and  we  get: 

fff 

A  =  a  -  =  1 
y)  b 


a  = 


Next,  look  at  the  variance: 

3  B  Qx-l02f(x)dx 
A  /“oofWdx 

Calculate  the  numerator  B  by  evaluating  the  integral 
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p  CO  n  CO 

B=  I  (x  —  p)2f(x)dx  =  I  a(x  —  [i)2e_b(-x_^2  dx 

j  —  CO  J  —  CO 

Use  the  change-of-variable  z=x-|j,  with  dz=dx,  x=-co— >z=-co,  and  x=co— »z=co,  to  get: 


B 


-■£ 


z  e 


dz 


Evaluate  this  integral  using  integration  by  parts: 


"  r 

udv  =  uv  — 


vdu 


with: 

u  =  z  -»  du  =  dz 

dv  =  ze_bz2dz  =  ( — — )  (— 2bz)e_bz2dz  -»  v  =  f — — ^ 
\  2b/  V  2b, 

The  integral  for  B  now  becomes: 


e 


-bz^ 


B 


=■£ 


z2e  bz2  dz 


=  (- 


3.Z  \  ,  2 

— ) e_bz 
2b^ 


+ 


2b 


£ 


e  bz2  dz 


The  first  part  of  the  integral  seems  to  be  the  product  (x, )(()),  but  we  can  eliminate  this  problem  by  using 
L’Hopital’s  Rule  as  follows: 


lim  ze 

Z— >oo 


d  f  \ 

z  Hi  ^  1 

lim  — ; — 7  =  lim  — ; - =  lim 


Z^CXD  gbz2  z— >00  d  /  ^z2  A  Z— >oo  2bzebz2  00 

di  > 


=  —  =  0 


The  same  thing  happens  when  z— >-oo,  so  the  first  part  of  the  integral  is  zero.  Except  for  the  constants  in 
front  of  the  integral,  the  second  part  of  the  integral  for  B  is  exactly  the  same  as  the  integral  for  A.  So, 
substituting  the  constant  a/2b  for  a,  we  get: 


B  = 


a  n 


2bA|b 

Using  the  values  for  A  and  B  in  the  equation  for  the  variance,  we  then  have: 
a  Fii 

„  B  2b  Jb  1 


G  A 


|rr  2b 
b 


b  = 


2a2 


Using  the  equation  for  a,  we  get: 

1  1 


b  _ 


tt  ^  2ttg2  V2ttg2 

So,  for  a  corresponding  normal  PDF,  the  values  of  the  coefficients  are  shown  in  Table  1: 
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Table  1:  Equations  for  the  Coefficients  of  a  Corresponding  Normal  PDF 


Coefficient 

Equation  for  the  Coefficient 

1 

a 

V2rr  ct2 

b 

1 

2o2 

c 

P 

Determining  the  Corresponding  Lognormal  PDF 

To  find  the  corresponding  lognormal  PDF  that  has  the  same  mean,  p,  and  variance,  a2,  of  a  given 
PDF,  we  can  start  by  writing  a  lognormal  equation  in  the  general  form: 

y  =  f(x)  =  ae_Kln(c))  (2) 

where  c  and  a  are  the  abscissa  and  ordinate,  respectively,  of  the  mode  of  the  corresponding  lognormal 
PDF  and  then  determine  the  value  of  b  so  that  the  area  under  the  normal  function  is  equal  to  one.  To  find 
the  values  of  a,  b  and  c.  Start  by  using  the  standard  equation  from  Calculus  for  finding  the  mean  x- value: 

N  J”xf(x)dx 
^  D  J”  f(x)dx 

D,  the  integral  in  the  denominator,  is  simply  the  area  under  the  lognormal  equation.  Calculate  the 
denominator  by  evaluating  the  integral: 


First  use  the  transformation: 


c 

x  =  cez 
with: 

dx  =  cezdz 

and  where  x=0— »z=-co  and  x=oo— >z=co.  The  transformed  integral  is: 

r°°  ,  ,  r°°  ,  , 

D  =  I  ae  bz  cezdz  =  ac  I  e  bz  +z  dz 

j  —  CO  J  —  CO 


Change  the  exponent  by  completing-the-square  to  get: 

-bz2  +  z  =  -b  (z2  -  -)  =  -b  (z2  -  -  +  — j - 2)  = 

V  b'  V  b  4b2  4b2/ 

-bz2  +  z  =  -b(z-— )  +  — 

V  2b/  4b 


so,  the  integral  now  becomes : 
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D  =  ac[ 


e-K-i) 

J-, 


dz 


Now  use  the  transformation: 
1 


m  =  z 


2b 


with  dm=dz,  and  where  z=-<x>— >m=-oo  and  z=co— >m=oo.  The  transformed  integral  is: 


D  =  ace4b 


i  r 
l-b 

« _ 


—  nnoih  I  zx  bm 


dm 


Except  for  the  constants  in  front  of  the  integral,  this  integral  for  D  is  exactly  the  same  as  the  integral  for  A 
in  the  previous  section  on  determining  the  corresponding  normal  PDF.  So,  substituting  the  constant 
aee1/l4b)  for  a,  we  get: 

D  =  ace4b 

If  the  corresponding  lognormal  curve  represents  a  normalized  PDF,  then  D=l: 

J_  /tc 

D  =  1  =  ace4b  — 

\  b 


Calculate  N,  the  numerator  of  the  mean  value  equation,  by  evaluating  the  integral: 
M  f  °°  -bfln(p))2 

N  =  axe  v  vc ))  dx 

40 

Again,  use  the  transformation: 


c 

x  =  cez 
with: 

dx  =  cezdz 

and  where  x=0— >z=-co  and  x=oo— >z=co.  The  transformed  integral  is: 

f  CO  f  CO 

N=  I  aceze_bz2  cezdz  =  ac2  I  e_bz2+2zdz 

j  —  OO  J  —  oo 


Change  the  exponent  by  completing-the-square  to  get: 


— bz2  +  2z  =  —  b  z 


— bz2  +  2z  =  — b  (z--)2  +- 
V  b/  b 


2z  1 

T  +  P 


The  integral  now  becomes : 
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N 


=  ac2  r 

J  —  c 


e  bvz  bJ  +b  dz  =  ac2eb 


1  r 
b 

j  —  c 


rb(z_&)2 


dz 


Now  use  the  transformation: 
1 

m  =  z  - 


b 

with  dm=dz,  and  where  z=-co— *m=-oo  and  z=co— »m=co.  The  transformed  integral  is: 


N  =  ac2eb 


£ 


e  bmZ  dm 


As  with  A  and  D,  again  use  the  standard  trick  to  evaluate  this  integral  and  get: 
7  JL  (if 

N  =  acre4b  — 

y  b 


So,  using  equations  N  and  D  and  the  equation  for  the  overall  mean  value,  jn,  we  get: 

ac2e4b 


A-  /tc 


N 

^  =  D 


ace4b 


A-  /tc 


b  A 

—  =  ce4b 


Next,  use  the  standard  equation  from  Calculus  for  finding  the  variance: 

2  _ M _ C(X  -  h)2f(x)dx 

D  J”  f(x)dx 

D,  the  integral  in  the  denominator,  is  simply  the  same  integral  as  we  evaluated  above  in  the  equation  for 
the  mean.  Calculate  the  numerator  M  of  the  equation  for  the  variance  by  evaluating  the  integral: 


M 


=  [  (x  —  p)2f(x)dx  —  [  (x  —  p)2a  e  b(ln©  dx 
Jn  J  n 


M 
with: 

Mi 


=  ■/" 
40 


x^  e 


-b 


dx  —  2ap 


f°xe-bK)/dx  + 

j  0  Jo 


dx 


=-/" 

40 


M2  =  2ap 


M3  =  ap‘ 


f™xe  b(ln©) 

40 

f”e-b('"©)2 

40 


dx 


Except  for  the  constants  in  front  of  the  integrals,  the  equations  for  M2  and  M3  are  exactly  the  same, 
respectively,  as  the  equations  N  and  D,  above.  Adjusting  for  these  different  constants,  we  get: 


M2  =  2pN  =  2pac2eb 
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i  Fa 

Mo  =  u2D  =  Li2ace4b  — 

\b 


To  solve  for  Mi,  we  again  use  the  transformation: 
c> 


=  '»© 


-  =  ez 
c 

x  =  cez 
with: 

dx  =  cezdz 

and  where  x=0— >z=-co  and  x=oo— >z=co.  The  transformed  integral  is: 


ra 

Mx  = 

J-  c 


ac2e2ze  bz2  cezdz  =  ac3 


ra 
J  —r 


-bz2+3z 


dz 


Change  the  exponent  by  completing-the-square  to  get: 


-bz2  +  3z  =  -b  (z2  -  v)  =  -b  (z2  - 1  +  ^  =  -b  ((2  - - 


4b^ 


3  \  ^ 

— bz2  +  3z  =  — b  (  z - )  + 


2b/ 

The  integral  now  becomes: 


Mj  =  ac3  j 


4b 


4b  dz  =  ac3e4b 


'£ 


b-KHe)2 


dz 


Now  use  the  transformation: 
3 


m  =  z  - 


2b 


with  dm=dz,  and  where  z=-<x>— >m=-oo  and  z=co— »m=oo.  The  transformed  integral  is: 


Mi  =  ac3e4b 


ra 

J  —  c 


—  I  A-bm2 


dm 


As  with  A  and  D,  again  use  the  standard  trick  to  evaluate  this  integral  and  get: 


,  JL  Iff 
Mi  =  ac  e4b  — 


Bringing  together  the  equations  for  Mi,  M2  and  M3,  we  get: 

,2L/fr  _  i  Iff  J_  Iff 

M  =  Mi  —  M2  +  M3  =  ac3e4b  j—  —  2pac2ebJ  — +  pzace4bj— 


b 


b 


2L  frf/  ,1  \ 

M  =  ace4b  I— I  c  eb  —  2pace4b  +  I 

The  coefficient  term  in  front  of  the  parentheses  is  equal  to  D.  Using  this,  we  get: 
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M  =  ace4b 

M  2  _3_ 

—  =  a2  =  c2eb  —  2[ice4b  +  q2 

Using  the  equation  we  derived,  above,  for  the  overall  mean  value,  q: 

3_ 

q  =  ce4b 


c2eb 


2qace4b  +  q' 


:)  =  D( 


c2eb 


2qce4b  +  q2 


we  get: 

g2  =  c2eb  —  2qce4b  +  q^  =  c^e 
1 


2  A  3  +  1  /  _3_\  J_  /  J_  \ 

AK  -1-  "2  -  -2"  2b  -  2[i(  ce4b  j  +  |12  =  |i2e2b  -  2 [i2  +  q2  =  q2  (e2b  -  2  +  1  j 


a2  =  [l2  (e2b  -  lj 

We  now  have  three  equations — for  |j,  o2  and  D=1 — to  solve  for  three  unknowns,  a,  b  and  c.  Starting  with 
the  equation  for  a2,  we  can  solve  for  b: 

g2  =  [i2  (e2b  -  lj 


a2  A 

—  =  e2b  -  l 
|i2 

A  a2 

e2b  =  —  +  l 

b2 


2b 
b  = 


=  In  — +  1 


21n 


Using  equation  for  |j,  we  get: 


q  =  ce4b 


c  =  qe  4b  =  qe  4' 


HSU)  =  ve'"&iV  =  J°± 


+  1  =n 


V-  +  °  '  (p2  +  a2)2 


The  equation  for  D=1  allows  us  to  calculate  the  abscissa  of  the  mode  of  the  corresponding  lognormal 
cost -risk  function,  given  the  function’s  mean  and  the  variance.  Using  the  equations  for  b  and  c,  we  get: 


A  fff 

D  =  1  =  ace4b  — 
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(p2  +  a2)2 

a  = _ p^ _  H2  +  a2 

So,  for  a  corresponding  lognormal  PDF,  the  values  of  the  coefficients  are  shown  in  Table  2: 

Table  2:  Equations  for  the  Coefficients  of  a  Corresponding  Lognormal  PDF 


Coefficient 

Equation  for  the  Coefficient 

p2  +  CT2 

a 

p3 

A 

l 

b 

2I„(?+1) 

F4 

c 

(p2  +  CT2)2 
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Addendum  2:  Mean  &  Variance  of  the  Convolution  of  Functions  of  Uncorrelated  Random 

Variables 


Determining  the  Mean  of  the  Convolution  of  Functions  of  Uncorrelated  Random  Variables 

If  all  the  functions  in  a  convolution  are  functions  of  uncorrelated  random  variables,  we  can  derive  an 
exact  relationship  between  the  mean  value  of  a  convolution  and  the  mean  values  of  the  individual 
functions  being  convolved.  Let  f(x)  be  the  convolution  of  two  functions  a(x)  and  b(x),  i.e., 
f(x)=a(x)*b(x),  and  let  these  two  functions  have  mean  values  |ja  and  |ib,  i.e.: 

fZ  xa(x)dx 
^  L”coa(x)dx 

fZ  xbQQdx 

^  J“00b(x)dx 

Since  the  integrals  in  the  denominators  are  simply  the  areas  A  under  the  functions  a(x)  and  b(x),  we  have: 


Ab  =  I  b(x)dx 

J  —  CO 


giving  us: 

J_0000xa(x)dx 


ha 


Aa 


xa(x)dx 


haAa  =  f 

j  —  CO 

and: 

fZ  xb(x)dx 


hb 


hbAb 


-L 


Ab 

CO 

xb(x)dx 


The  mean  value  |if  of  the  convolution  f(x)  is: 

fZ  xf(x)dx 


hf  = 


IZ  fWdx 


Since  f(x)  is  the  convolution  of  a(x)  and  b(x),  then  the  area  under  f(x)  is  the  product  of  the  areas  under 
a(x)  and  b(x),  i.e.,  Af=AaAb,  so  we  then  have: 


bf : 


fZ  xf(x)dx 
Af 

xf(x)dx 


hfAf  =  J 

and: 
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p  LXJ  p  CU  p  UU  pU U  /iUO 

UfAf  =  xf(x)dx  =  x  a(x  —  a)b(a)dadx  =  xa(x  —  a)b(a)  dxda 

J  —  CO  J  —  CO  J  —  CO  J  —CO  J  —  00 


By  using  p=x-a  (— >x=a+P)  and,  if  we  take  the  derivative  of  p  with  respect  to  x  and  express  the  result  in 
differential  form,  we  get  dp=dx.  We  then  have: 

p  CO  p  CO  p  CO  p  CO 

PfAf  =1  I  (a  +  (3)a((3)b(a)  d(3da  =  I  I  (aa(p)b(a)  +  (3a((3)b(a))  d(3da 

j  —  CO  J  —  CO  J  —  CO  J  —  CO 


PfAf  - 


aa((3)b(a)  d(3da  + 


(3a((3)b(a)  d(3da 


If  the  variables  a  and  p  are  uncorrelated  random  variables,  then  the  functions  a(x)  and  b(x)  are  separable, 
giving  us: 


p  CO  /  p  CO  \  p  CO  /  p  CO  \ 

PfAf  =  I  a((3)  (  I  ab(a)dajd(3  +  I  b(a)  |  I  pa(p)dp|da 

j  —  CO  \j —CO  /  J  —  CO  \J —  CO  / 

p  CO  /  ^  CO  \  p  CO  /  ^  CO 

PfAf  =  I  a(p)dp(  I  ab(a)daj+  I  b(a)da|  I  (3a(p)d(3 

j  —  CO  \J  —  CO  I  J  —  CO  \  J  —  CO 


The  quantities  in  the  parentheses  are  the  mean  values,  and  the  quantities  in  front  of  the  parentheses  are  the 
areas  Aa  and  Ab.  So,  the  terms  being  added  are  mean  values  and  their  associated  areas,  paAa  and  pbAb- 
giving  us: 

PfAf  —  paAaAb  +  pbAaAb  —  (pa  +  pb)AaAb 


But,  since  Af=AaAb,  and  we  have: 


hf  =  ha  +  hb 


So,  for  functions  of  uncorrelated  random  variables,  the  mean  of  a  convolution  is  the  sum  of  the  means  of 
each  of  the  functions  in  the  convolution.  If  the  functions  a(x)  and  b(x)  are  themselves  convolutions  of 
functions  of  uncorrelated  random  variables,  then  we  can  express  each  of  their  means  as  the  sum  of  the 
means  of  the  functions  in  their  convolutions  and  we  can  then  express  the  mean  of  f(x)  as  the  sum  of  both 
sets  of  means.  We  can  extend  this  indefinitely  for  any  set  of  n  uncorrelated  functions  to  get: 


Pf  =  Pi  +  P2  +  Ps  +  ...+  pn 


I" 


where  n  is  any  positive  integer. 


Determining  the  Variance  of  the  Convolution  of  Functions  of  Uncorrelated  Random 
Variables 

For  functions  of  uncorrelated  random  variables,  we  can  also  derive  an  exact  relationship  between  the 
variance  of  a  convolution  and  the  variances  of  the  individual  functions  being  convolved.  Let  f(x)  be  the 
convolution  of  two  functions  a(x)  and  b(x),  i.e.,  f(x)=a(x)*b(x),  and  let  these  functions  have  mean  values 
o2f,  o2a  and  o2b,  respectively.  The  variance  o2f  of  the  function  f(x)  is: 


at  = 


J”  (x  -  pf)2f(x)dx  /“  (x  -  pf)2f(x)dx 


Lo 0  fWdx  Af 

p  CO  p  CO  p  CO 

AfG2  =  I  (x  —  pf)2f(x)dx  =  I  (x— pf)2  I  a(x  —  a)b(a)dadx 

J —  CO  J  —  CO  J —  CO 
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Affff  = 


(x  —  |if)2a(x  —  a)b(a)  dxda 


By  using  p=x-a  ( — >x=</+[i)  and,  if  we  take  the  derivative  of  p  with  respect  to  x  and  express  the  result  in 
differential  form,  we  get  dp=dx.  We  then  have: 

n  CO  CO 

Af(Jf  =|  I  (a  +  p  -  pa  -  pb)2a(p)b(a)  dpda 

j  —  CO  J  —  CO 
n  CO  ^  CO 

Affff  =|  I  (a2  +  p2  +  pa2  +  pb2  +  2«p  -  2apa  -  2apb  -  2Ppa  -  2Ppb  +  2papb)a(p)b(a)  dpda 
j  —  00  *2  —  00 

^  CO  ^  CO 

Affff  =  I  I  ((P2  -  2Ppa  +  pa2)  +  (a2  -  2apb  +  pb2) 

*2  —  00  J  —  CO 

+  2(«p  -  apa  -  Ppb  +  papb))a(p)b(a)  dpda 

^  CO  ^  CO 

Affff  =|  I  (P  —  pa)2a(P)b(a)  dpda 

J  —  CO  *2  —  00 

p  CO  p  CO 

1  I  (a  -  pb)2a(p)b(a)  dpda 

j  —  CO  J  —  CO 
n  CO  p  CO 

2  I  I  (a p  -  apa  -  Ppb  +  papb)a(p)b(a)  dpda 

*2—00  J  —  CO 
n  CO  r-  CO 

Affff  =1  (P  -  pa)2a(p)dp  I  b(a)da 

*2  —  CO  *2  —  CO 

/*  CO  /-*  CO  n  CO  n  CO 

+  I  (a  —  pb)2b(a)da  I  a(p)dp  +  2  I  ab(a)da  I  pa(p)dp 

*2  —  CO  *2  —  00  *2  —  CO  *2  —  00 

p  CO  CO  /*  CO  /i  CO 

+  2pa  I  ab(a)da  I  a(p)dp  +  2pb  I  b(a)da  I  pa(p)dp 

J  —  CO  *2  —  CO  *2  —  CO  J —  CO 

nCO  n  CO 

+  2papb  I  b(a)da  I  a(p)dp 

*2  —  CO  *2  —  CO 


+ 

+ 


Affff  —  AaffaAb  +  AbffbAa  +  2(AapaPbAb  paPbAbAa  pbpaAaAb  +  papbAaAb) 

Afff2  =  fflAaAb  +  ffbAaAb  +  2(0)  =  ffiAaAb  +  ffbAaAb 

Since  f(x)  is  the  convolution  of  a(x)  and  b(x),  then  the  area  under  f(x)  is  the  product  of  the  areas  under 
a(x)  and  b(x),  i.e.,  Af=AaAb,  so  we  then  have: 


AaAbfff  —  ffaAaAb  +  ffbAaAb 

fff  —  aa  +  ab 


So,  for  functions  of  uncorrelated  random  variables,  the  variance  of  a  convolution  is  the  sum  of  the 
variances  of  each  of  the  functions  in  the  convolution.  If  the  functions  a(x)  and  b(x)  are  themselves 
convolutions  of  functions  of  uncorrelated  random  variables,  then  we  can  express  each  of  their  variances 
as  the  sum  of  the  variances  of  the  functions  in  their  convolutions  and  we  can  then  express  the  variance  of 
f(x)  as  the  sum  of  both  sets  of  variances.  We  can  extend  this  indefinitely  for  any  set  of  n  uncorrelated 
functions  to  get: 


ff2f  =  ff2l  +  ff22  +  ff2  3  +  ...  +  ff2n 


n 


i=l 


where  n  is  any  positive  integer. 
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Determining  the  Mean  and  Variance  of  the  Convolution  of  Discrete  Functions 

In  the  special  cases  of  discrete  functions,  where  the  functions  are  a  set  of  one  or  more  delta  functions,  we 
can  use  the  convolution  properties  of  delta  functions  to  see  how  they  affect  the  mean  and  the  variance.  If 
f(x)  is  the  convolution  of  any  real  function  a(x)  with  a  delta  function  b(x)=8(x-x0),  we  get: 

pOO  p  CO 

f(x)  =  a(x)  *  b(x)  =  I  a(a)b(x  —  a)da  =  I  a(a)5(x  —  a  —  x0)da 

j  —  CO  J  —  CO 

Let  p=x-x0: 

p  CO  p  CO 

f(x)  =  I  a(a)5(x  —  a  —  x0)da  =  I  a(a)5((3  —  a)da  =  a(p)  =  a(x  —  x0) 

j  —  CO  J  —  CO 

So,  a(x)  is  simply  shifted  to  the  right  by  the  amount  x0.  This  will  also  shift  the  mean  to  the  right  by  the 
amount  x0  and  leave  the  variance  of  a(x)  unchanged.  So,  if  f(x)  is  the  convolution  of  any  real  function 
a(x)  with  a  set  of  n  delta  functions: 


11 

b(x)  =  ^  q  5(x  -  x;) 

i=l 


where  n  is  any  positive  integer,  the  Cj  are  constant  coefficients  and  the  Xi  are  fixed  real  numbers.  We  then 
get: 


f(x)  =  a(x)  *  b(x)  =  j  a(a)b(x  —  a)da  =  J  a(a)  q  5(x  —  x;) 

=  ^qj  a(a)5(x  —  a  —  X;)da 
i=l  -co 


I  da 


f(x) 


Let  Pi=x-Xj: 


f(x) 


Zp  CO  _  P  CO  _  _ 

q  I  a(a)5(x  —  a  —  X;)da  =  y  q  I  a(a)5((3i  —  a)da  =  y  q  a(P;)  =  y  qaCx  —  Xj 

„  J  —  CO  (  J  J  —  CO  .  *  (  !* 

1  =  1  1  =  1  1  =  1  1=1 


i) 


So,  f(x)  is  the  sum  of  n  copies  of  a(x)  with  each  ith  copy  shifted  to  the  right  by  the  amount  Xj.  The  mean  p 
of  f(x)  is: 


h 


J_°L  xf(x)dx  _  x  SjLi  q  a(x  -  Xj)dx  _  Ef=1  q  xa(x  -  xt)dx 
OMdx  Ef=1  C;  a(x  —  Xj)dx  Sf=1q  J_0000  a(x-xi)dx 


The  denominator  is  the  area  Af  under  the  function  f(x).  Af  is  the  weighted  sum  of  the  areas  Aa  under  the  n 
shifted  copies  of  the  function  a(x): 


ii  ii 

Af  —  q  Aa  —  Aa  C; 

i=l  i=l 

Using  this,  we  get: 


Afp 


n  n  oo  n  n 

=  Aa  ^  q  =  ^  q  J  xa(x  -  x;)dx  =  ^  q  Aap;  =  Aa  ^  q  p .j 

i= 1  j— ^  — 00  y=i  i=l 
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Aa  V 


So  the  mean  ju  of  f(x)  is  the  weighted  average  of  the  means  of  the  n  shifted  copies  of  a(x). 


The  variance  o'  of  f(x)  is: 

J_°L(x  -  p)2f(x)dx  J_>  -  p)2(Ef=i  q  a(x  -  xj))dx  Z"=1  q  J_°^(x  -  p)2a(x  -  x;)  dx 


9  J  —  00 

V“  -  - 


f™oo  f(x)dx 


Af 


Af 


n 


2  _Xf=1CiAaG2i  _  Aav- 
a  Af  Ci  °  1 


Af 


i=l 


So  the  variance  of  f(x)  is  the  weighted  average  of  the  variances  of  the  n  shifted  copies  of  a(x). 
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